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ABSTRACT 

The Abelian Higgs model with anisotropic couplings in 2 + 1 dimensions is 
studied in both the compact and non-compact formulations. Decoupling of the 
space-like planes takes place in the extreme anisotropic limit, so charged particles 
and gauge fields are presumably localized within these planes. The behaviour of 
the model under the influence of an external magnetic field is examined in the 
compact case and yields further characterization of the phases. 



1 Introduction 



The existence of a layered phase in gauge theories with anisotropic coupling has 
been first conjectured by Fu and Nielsen ([|Ij], see also 0) in the early eighties. 
If the couplings are equal in all directions, the Wilson loops are governed by the 
area law (if the couplings are strong) or by the perimeter law (if they are weak). 
If now the couplings are weak in d dimensions and strong in the remaining one 
dimension, then the model will exhibit Coulomb forces in the d-dimensional sub- 
space and will confine in the last direction; this means that charged particles will 
be localized within the d-dimensional subspace. In the extreme case, where the 
coupling in the extra dimension is infinite, charged particles from each subspace 
will remain strictly confined to the d-dimensional space-time; in other words, 
the d-dimensional hyperplanes(layers) will decouple from each other and this is 
at the origin of the characterization of this phase as "layered" . 

Monte Carlo evidence for a layered structure in pure U(l) has been given 
later on |3j . There exists a lowest dimension for which the layered phase exists: 
it is the dimension above which the system has two phases. For the abelian 
gauge theory it is at five dimensions that one may first discriminate between 
the Confinement and the Coulomb phases; for the non-abelian case the critical 
dimension is six. 

The initial motivation for the study of a model with anisotropic couplings has 
been dimensional reduction with cosmological implications: if someone lives in 
one of the d-dimensional hyperplanes of a system in the layered phase, he has 
no way of communicating with the remaining layers and a kind of dimensional 
reduction is achieved. These considerations have also been given a modern version 
in the context of membrane theories and confining (or localization) of particles 
within the branes f|. On the other hand, the construction of the chiral layer is 
an important ingredient of the models of chiral fermions on the lattice ||. 

In this paper we try to perform an exploratory study of a richer model. We 
treat a (2+l)-dimensional Abelian Higgs model. This has the advantage that it 
is a low- dimensional model with a phase transition (at large gauge couplings), 
so it may have a layered phase. In addition, this system has a direct relation 
to the ceramics exhibiting high-T c superconductivity . These materials have a 
layered structure, so they are effectively planar and high-T c superconductivity is 
considered as a two dimensional phenomenon. In this sense the (2 + 1) - dimen- 
sional theory acquires physical relevance. The isotropic version of the model has 
been extensively studied recently P, ||, in connection with cosmological consid- 
erations and Condensed Matter Physics, in which external magnetic fields play a 
predominant role, such as penetration and vortex formation. 

Another feature of the various phase transitions associated with the model is 
that they may be strengthened in the presence of an external magnetic field; in 
addition, the phases of the model with compact action may be characterized by 
the behaviour of the system under the influence of the magnetic field. 

Both the compact and non-compact versions of the model have been examined 



2 



in this paper and they both have characteristics that one can attribute to layer 
decoupling. 



2 Formulation of the model 

The model under study is the Abelian Higgs model in the three-dimensional space. 
Direction 3, corresponding to z, will be singled out by couplings that will differ 
from the rest; later on an external magnetic field will also point in this direction. 

We first discuss the non-compact version of the model; we proceed with writing 
down its lattice action. 

S = ~J gs E + + FUx)} 

Z X Z X 

+\(3hs (^2[2(p*(x)ip(x) - y*{x)Ui{x)y{x + i) - ip*(x)U^x)ip(x + 2)] + h.c^j 

(Tlv'ixMx) - <P*(xM(xMx + 3)] + h.c^j 
+ £[(1 " Wr - Whs - p M )<p*(x)<p{x) + p R (tp*(x)tp(x)) 2 } (1) 

X 

where F^- = Aj(x + i) — Aj(x) — Ai(x + j) + Aj(x). 

We have allowed for different couplings in the various directions: the ones 
pertaining to directions 1, 2 are given the subscript "s" for "space", while direc- 
tion 3 gets couplings with the subscript "t" for "time" . We stress that the third 
direction is not really "time"; it is just the direction singled out as explained 
previously. However, we use this notation for convenience; in particular, we refer 
below to "time-like" plaquettes, links or couplings with this understanding. 

For future use we mention here some notations. The link variables U^(x) 
are defined as e %9sCXsAs or e %9tatAt respectively. They are also written in the form 
U^{x) = e l6 ^ x \ since they are complex phases. In addition, the scalar fields are 
also written in the polar form (f(x) = R(x)e lx<yX \ The order parameters that we 
will use are the following: 

Space - like Plaquette : P s =< ^ F n( x ) > ( 2 ) 

^ ' X 

Time - like Plaquette : p t =<^ Y.( F U*) + ^0*0) > (3) 
Space - like Link : L s =< — L ]T[ C os(x(:r + i) + B^x) - X (x)) (4) 



X 



+ cos(x(a: + 2) + 6$(x) - x(x))} > (5) 
Time - like Link : L t =< £ cos( X (x + 3) + 6${x) - x(x)) > (6) 
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Higgs field measure squared : p 2 



N 3 



(7) 



We now proceed to the naive continuum limit (at tree level) of the lattice 
action [10]. The first step is to rewrite the action in terms of the continuum fields 
(denoted by a bar): 



^) 1/2 , 

Phs 



Ax = a 8 A x , A 2 = a s A 2 , A 3 = a t A 3 . 

This means that the time-like part of the pure gauge action is rewritten in the 
form: f3 gt a t J2a 2 s a t [(F 13 ) 2 + (F 23 ) 2 } -> f3 gt a t J d 3 x[(F 13 ) 2 + (F 23 ) 2 }. On the other 

hand the space-like part is: (3 gs ^I2 a 1 a t(Fi 2 ) 2 



P 9S ^Jd 3 x(F 12 ) 2 . If we define 



at 

dWs' 



P 9 t = 



gfat ' 



the resulting continuum action reads 

1 



d 6 x 



9t 



(F 12 ) 2 + -[(F 13 y + (F 



9t 



2:i) 



gs_Os 



Defining 7 5 = {^Y^ 2 and using the definitions of /3 gs , (3 gt we find that 7 9 
We denote by £ the important ratio ^ of the two lattice spacings (the correlation 
anisotropy parameter) and finally derive the relation: 



(1/2) 



9i 



Along the same lines, one may rewrite the scalar sector of the action in the form: 



J d 3 x[\D^\ 2 + \D 
We have used the notations: 



2<P\ 2 



7, 



^\D 3 <p\ 2 + m 2 rv + HrW] 



Jhs, 



1/2 



2 2 

m a. 



il-2[3 R -2(3 hs -p ht ), \a s 



4(3 



R 



Phs Phs** 

Of course, the dependence of the space-like and time-like couplings on the 
lattice spacings a s and a t is a very interesting issue, which has been pursued 
in the context of QCD flO| , 11 , 12| . Neglecting quantum corrections we have 
1^ = l g = £• In this paper we would like to explore the decoupling of the space- 
like planes, so in principle we should have a relation for the variation of £ as a 
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function of 7^ and 7 g , with the quantum corrections properly taken into account. 
We defer this study to a later stage and for the time being we adopt the tree 
level scaling refered to previously. In particular, we will vary the quantities 
/3 gs , Pgt, Phs, Pht in such a way that we have 

l v = 7 9 = C- 

We have used the letter ( for the common ratio of the couplings, so our assumption 
reads: f = ( and g s = g t = g. 

The parameter /3r is found from the equation (3r = using the fixed value 

x = 2 for the parameter x=\. 

The compact action is defined as: 



Scompact = figs XX 1 ~~ COsF 12 (x)) + (3 gt K 1 ~ COsF 13 (x)) + (1 - COsF 23 (x))] 

x x 

+f3 hs Yftv\x)v(x)-v\x)Ui(x)e iA ^ {x) ^ 

X 

+/3 ht J2&*(xMx) - ip*{x)U % {x)e iA ^W<p{x + 3)] 

X 

+ - 2(3 R - 2(3 hs - /3 M )(p*(x)<p(x) + p R (<f*{x)<f(x)) 2 } (9) 

X 

The difference from the non-compact version lies in the gauge kinetic terms; 
however, since we will treat the compact system in an external magnetic field, we 
have taken the opportunity to also include in the action the background potential 
^■kexti x )i which will generate the external field. There are many ways to impose 
such an external magnetic field [|X3j, |], [TJ]]; for example, one might include it only 
in the interaction term with the matter fields; that is, the matter fields interact 
with both the quantum and the background fields. This has been our choice, as 
can be seen from equation (|9|). 

The order parameters in this case are defined in a very similar way as in the 
non-compact case. We state the expressions which differ: 

Space — like Plaquette : P s =< — - cosF 12 (x) > (10) 

^ * X 

Time - like Plaquette : P t =< — — ^ ^2[cosF 13 (x) + cosF 23 (x)] > (11) 



3 Algorithms 

We used the Metropolis algorithm for the updating of both the gauge and the 
Higgs field. It is known that the scalar fields have much longer autocorrelation 
times than the gauge fields. Thus, special care must be taken to increase the 



5 



efficiency of the updating for the Higgs field. We made the following additions to 
the Metropolis updating procedure ||: 

a) Global radial update: We update the radial part of the Higgs field 
by multiplying it by the same factor at all sites: R(x) — > e^R(x), where £ G 
[— e, e] is randomly chosen. The quantity e is adjusted such that the acceptance 
rate is kept between 0.6 and 0.7. The probability for the updating is P(£) = 
min{l, exp(2V£ — A£>(£))} where AS(£) is the change in action, while the 2V^ 
term comes from the change in the measure. 

b) Higgs field overrelaxation: We write the Higgs potential at x in the 
form: 



V(<p(x)) = -a F + R 2 (x) + (3 R (R 2 (x) - 1) 



(12) 



where 



a = 



R(x) cosx(^) 
R(x) sinx(x) 



F x = f3, 



h a 



[Rip + «) cos(x(x + + 6i(x)) 



i=l,2 



P, 



1 1 6' 



^2 R( x + i) sin(x(x + i) + 0i(x)) 



i=l,2 



+ p ht R(x + 3) cos(x(x + 3) + 9$[x) ) , 



+ p ht R{x + 3) sin(x(x + 3) + 6$(x)). 



We can perform the change of variables: (a, F) — > (X, F, Y) , where 



F = IFI 



X = a ■ f , Y = a — Xf . 



(13) 



The potential may be rewritten in terms of the new variables: 

V(X, F, Y) = -XF + (1 + 2(3 R (Y 2 - 1))X 2 + Y 2 (l - 2f3 R ) + (3 R {X 4 + Y 4 ). (14) 

The updating of Y is done simply by the reflection: 

Y -f Y' = -Y. (15) 

The updating of X is performed by solving the equation: 

V(X',F,Y) = V(X,F,Y) (16) 

with respect to X'. Noting that X' = X is obviously a solution, we may factor 
out the quantity X' — X and reduce the quartic equation into a cubic one, which 
may be solved. The change X — > X' is accepted with probability: P(X') = 
min{P , 1}, where P - 9 -^£^/ dV %£X) . 
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4 Results 



4.1 Runs with fixed £ 

In this set of measurements we set f3 gs = 6.0 ,x = 2, and let Ph s run. The 
remaining coupling constants vary according to the value of (, that is: 

Pht = ( 2 Phs, P g t = CP 9 s, PH=^T- (17) 



4/9, 



We recall that C is the ratio: C = \ r = \ fn - in figure II] we show a comparison 

y ^>ff s V Phs — 

of the behaviour of the measure of the scalar field p 2 for the symmetric (( = 1.0) 
and a highly asymmetric model (( = 0.1), that we call "S(ymmetric) model" 
and "A(symmetric) model" respectively. We have used the compact action in 
both models. The S model (( = 1.0) shows that the system undergoes a phase 
transition starting at Phs — 0.35 from the 3D Coulomb phase to the 3D Higgs 
phase. (We note here that p 2 is not going to infinity for big Ph s , which is the 
case in several treatments; the reason is that /3r is not fixed, but increases with 
increasing /3 hs ). Then there is a transition region from p hs ~ 0.35 to Ph s — 0.50; 
for /3hs > 0.50 the system is in the Higgs phase, characterized by a large value of 

The corresponding transition for the A model (( = 0.1) moves towards bigger 
values of Ph s arid it seems smoother. In the latter case the "time-like" coupling 
constants are ( 2 = 0.01 times smaller than their "space-like" partners. Let us 
note that for ( = there will be no communication at all among the planes. 
This lends support to the assumption that for ( = 0.1, that is ( 2 = 0.01, we may 
have reached the limit of decoupled layers; indeed, we have checked that for an 
even smaller ratio of couplings (( = 0.01) the curve for the A model does not 
change much. This presumably means that the time-like separation a t of the 
spatial planes is already much bigger than the spatial lattice spacing a s . This 
is further supported by the measurement of other quantities, which depend on 
the direction (see, for example, figure 2 and especially figure 3 below). It turns 
out that the quantities P s , L s are much bigger than P t , L t for all the range of 
Phs, indicating that the quantities related to the communication of the planes are 
negligible in this compared against the similar quantities within the layers. 

It seems safe to assume that the region Phs < 0.55 corresponds to a Coulomb 
phase, where the layers are decoupled: this is equivalent to the Coulomb phase 
of the corresponding two-dimensional model. The transition region extends from 
P hs ~ 0.55 to Phs — 0.65. and then we have a Higgs phase for the A model: it 
has smaller p 2 than the one of the S model and the quantities related to the third 
dimension are very small. This presumably means that the layers are decoupled 
also in this phase, so the picture is that we have moved effectively to the Higgs 
phase of the corresponding two-dimensional model. Thus we may say in brief that 
in the A model we see a transition from the 2D Coulomb phase to the 2D Higgs 
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< R A 2 >, bg_sp=6.0, x=2.0, vol=8 A 3 



zeta=1 .0 
zeta=0.1 



* + X X X X X x 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 

betahsp 



Figure 1: p 2 versus flhs for the S model ( ( = 1.0) and the A model (£ = 0.1) for 
a 8 3 lattice with f3 gs = 6.0 and x — 2. 

phase. Let us stress here that we have been using the words Higgs "phase" and 
"phase" transition in a rather loose meaning when refering to the asymmetric 
model. We have not checked that the passage from one region to the other is 
actually a phase transition: it may be a crossover, so the exact characterization 
is open for the time being. We should remark that the critical parameter for the 
Higgs phase transition is of order |, where d is the space-time dimension. Thus, 
it is not accidental that the isotropic model has a phase transition at (3h s — 0.35; 
this is close to the expected value, since d=3. On the other hand, the anisptropic 
model is effectively two-dimensional, so one should expect a value approximately 
equal to |, in semi-quantitative agreement with our results. 

We note that the relative position of the transitions means that for f3 hs < 0.35 
the system lies in the Coulomb phase for either value of (; for j3h s = 0.50 the S 
model lies in its 3D Higgs phase, while the A model is still in its 2D Coulomb 
region; for (3}- ls = 0.70 both models lie in their respective Higgs phases. The region 
between the two curves of figure 1 is presumably full of other curves corresponding 
to all values of ( between 0.1 and 1. Later on we will consider the behaviour of 
several quantities with fixed flhs but varying (; figure 1 suggests that we may use 
fihs = 0.5 and j3h s = 0.7, based on the comments we made above. 

Figure 2 shows the same models as above, but the space-like link is depicted 
rather than p 2 ; in addition, we have shown the results for both the compact 
and the non-compact version of the model. The changes take place at the same 
values of f3 hS) as in figure 1. The results of the non-compact model are very close 
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Link sp, bg_sp=6.0, b_t=b_sp*zeta A 2, vol=8 A 3 



0.9 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



* 



x zeta=1 .0, non-compact 

+ B compact 

0.8 I- x v zeta=0.1 , non-compact 



compact □ 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

betahsp 



Figure 2: L s versus fihs for the S model (£ = 1.0) and the A model (( = 0.1) for 
a 8 3 lattice with f3 gs = 6.0 and x = 2.0. 

to the ones of the compact model, the most sizable differences being observed in 
model S, where the discrepancy reaches the value 0.2 in the region of the phase 
transition. In model A the differences are really small. 

In figure 3 we display the time-like link and plaquette only for the compact 
version of model A (£ = 0.1), to get a more concrete idea of the layered phase. 
The parameters vary in the same way as above. We find that the time-like link 
starts increasing around j3h s = 0.6, which is in the region where the passage in 
model A is located in figures 1 and 2. The time-like plaquette Pt does not seem 
to "feel" the passage and assumes a constant value. But the most important 
characteristic is that both quantities are really very small: their maximal value 
is 0.03, while their space-like partners L s and P s take on much larger values: 
0.13 < L s < 1.0, 0.90 < P s < 1.0. This is of course a consequence of our choice of 
bare parameters and confirms the fact that quantities related to communication 
between the planes are almost negligible in this region of the parameter space. 
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0.035 



0.03 



0.025 



0.02 



0.015 



0.01 
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bg_sp=6.0, Asymmetric (zeta=0. 10), vol=8 A 3, Compact U(1) 
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Link t x 



0.3 



0.4 
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betahsp 



0.8 



0.9 



Figure 3: Time-like link and plaquette versus f3h s for the A model (( = 0.1) for 
a 8 3 lattice with (3 gs = 6.0 and x = 2.0. 



It is possible to show in a very explicit way the decoupling of the planes if 
we restrict ourselves to the Coulomb phase. The Wilson loops W(N, N) can be 
calculated analytically in the 2D (pure) U(l) model and the resulting values can 
then be compared to the (space-like) plaquettes found by Monte Carlo simulations 
of the asymmetric Higgs model in the extremely layered Coulomb phase. In the 
table that follows we have put the results for the Wilson loops with sizes N x N 
from N = 1 to N = Q. The values in parentheses are the statistical errors in the 
last one or two digits. The values of the parameters are reported in the table 
and the quantity ( is set to 0.1. We also present in the table the values W(N, N) 
from the analytical expressions |0|: W(N, N) = (^jj) {NxN) - It is obvious that 
the system has moved to essentially two-dimensional dynamics, since most of the 
discrepancies lie within the error bars of the Monte Carlo data. 
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Table 1 


Volume = 


12 3 , C = 0.1, (3 hs = 0.30, x = 2.0 


Loop size 


Measured 


Analytic (2D, JT|) 






2.0 


1 x 1 


0.6981(1) 


0.697775 


2x2 


0.2380(2) 


0.237061 


3x3 


0.0395(2) 


0.0392136 


4x4 


0.0033(2) 


0.00315823 


5x5 


0.0002(2) 


0.000123845 


6x6 


0.0001(1) 


0.0000023645 




P 9 s = 


4.0 


1 x 1 


0.8635(1) 


0.863523 


2x2 


0.5562(2) 


0.556026 


3x3 


0.2669(4) 


0.266971 


4x4 


0.0955(5) 


0.0955827 


5x5 


0.0251(4) 


0.0255178 


6x6 


0.0053(2) 


0.00507988 




figs = 


6.0 


1 x 1 


0.9123(1) 


0.912359 


2x2 


0.6931(2) 


0.692889 


3x3 


0.4384(4) 


0.438019 


4x4 


0.2307(9) 


0.230491 


5x5 


0.1008(13) 


0.10096 


6x6 


0.0360(11) 


0.0368106 






8.0 


1 x 1 


0.9352(1) 


0.935235 


2x2 


0.7655(1) 


0.76504 


3x3 


0.5490(3) 


0.54738 


4x4 


0.3454(6) 


0.342559 


5x5 


0.1909(9) 


0.18751 


6x6 


0.0925(10) 


0.089775 



4.2 Runs with varying £ 

Now we come to a somewhat more detailed examination of the phase diagram 
where the "space-like" coupling constants are kept fixed and the ratio ( changes, 
also forcing (3 gt and (3ht to change. Of course, the result will strongly depend on 
the region in which the fixed "space-like" coupling constants lie; we will present 
results for the two interesting choices that we described in the discussion of figure 
1 ; Pgs = 6.0 and the values 0.5 and 0.7 for {3 hs . 

In figure f| we show the space-like link L s for the case with (3h s — 0.5; again 
both compact and non-compact results are displayed. The parameter ( starts 
from zero, where we expect a Coulomb phase with fully separated layers, that 
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Link sp, bg_sp=6.0, bh_sp=0.500, b_t=b_sp*zeta A 2, vol=8 A 3 



1 1 


1 1 1 1 1 

non-compact + 
compact x 




* * * ; 

. * * * 
* 

X 

X 


- 


X 

+ 

X 


- 


+ 

X 

+ 




X 




+ 




+ X 




+ + x 

j- + i x 
* X x x x x x 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

zeta 

Figure 4: L s versus ( for j3 gs = 6.0 and flhs = 0.500. 

is a 2D model, to 0.8, which corresponds to a model close to the symmetric 
3D one. The differences between the compact and the non-compact version are 
localized in the transition region. The non-compact results are somewhat bigger, 
the maximum difference being about 0.05, but qualitatively the results are similar 
for both versions. The transition from the symmetric to the broken phase takes 
place at about ( ~ 0.5, where 2D Coulomb becomes 3D Higgs: The behaviour 
of the space-like link may be interpreted with the help of figure 2. As we change 
(, we move continuously to new curves, of the sort of the ones in figure 2, but 
characterizing the £ under study. When ( = 0, the corresponding curve (close to 
the lowest one in figure 2) will inform us that the value flhs = 0.5 corresponds 
to the Coulomb phase for this value of (. The curves for larger ( will differ from 
the initial one in an important aspect: the transition region will move to the 
left (to smaller j3h s ) and will be steeper, until finally they will coincide with the 
uppermost curve in figure 2. When the transition region reaches f3h s = 0.5, the 
value of L s will increase abruptly and the system will pass over to a Higgs region. 
It is not clear at this point whether this is 2D or 3D Higgs, but this will be 
clarified by figure 5, which follows. 

The time-like link L t , shown in figure 5 both for the compact and non-compact 
versions, seems a more sensitive detector of the phase transition: it changes 
already in the region < ( < 0.5, where the space-like link is essentially constant, 
and its variation is bigger, since it starts from zero. Of course, it has to become 
equal to the space-like link (about one) in the limit ( — * 1, since there one 
regains the symmetric model. The non-compact results lie systematically above 
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LinkJ, bg_sp=6.0, bh_sp=0.500, b_t=b_sp*zeta A 2, vol=8 A 3 
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zeta 



Figure 5: L t versus ( for (3 gs = 6.0 and (3h s = 0.500. 



the compact ones. The comparison of this figure with the previous one shows that 
the time-like link also gets its 3D value at the same time when the space-like 
link jumps. This fact indicates that we have a transition from the 2D Coulomb 
phase to a 3D Higgs, that is, both changes (in dimension and in the breaking of 
symmetry) appear to happen simultaneously. 

In figure ^| we show the space-like link L s versus ( for the second of the two 
chosen values, (3h s = 0.7. In this case (5h s lies to the right of the transition already 
for the curve corresponding to ( = 0, as may be seen in figure 2. Thus we do 
not expect a phenomenon as in figure 4, where the transition region was passing 
through /3hs- We start from a layered Higgs phase already at ( = and we may 
isolate the transition from the 2D to the 3D Higgs phase (the space-like link has 
big values, from 0.90 to 0.97, in all the range of £, in striking contrast with figure 
4). The transition happens at f3h s ~ 0.25 for the compact case. It seems smooth, 
however one may observe a volume dependence in the results of figure 6, where 
we have put mesurements from both 8 3 and 16 3 lattices. In addition, we remark 
that the points for the compact 8 3 have a non-monotonous region between 0.15 
and 0.25; this disappears in the 16 3 lattice, so we interpret it as finite size artifact. 
The same effect is seen in the non-compact 8 3 case, but seems unimportant, in 
view of the experience with the compact case. The behaviour of the non-compact 
8 3 is smoother than its compact partner. 

Figure 7 contains the results for the time-like link L t versus ( for the two 
volumes 8 3 and 16 3 . This quantity is a good indicator of the 2D to 3D phase 
transition, as we said above; its volume dependence is almost negligible, but 
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Linksp bg_sp=6.0, bh_sp=0.700 b_t=b_sp*zeta A 2 
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Figure 6: Compact and non-compact L s versus ( for (3 gs = 6.0 and (3hs = 0.700. 

the difference between the compact and non-compact version is very pronounced 
at the region of the transition. The non-compact version appears steeper, in 
contrast with the previous figure 6. 

4.3 Vortex dynamics in the compact version 

We will now consider the system under the influence of a homogeneous external 
magnetic field; thus one should construct first a lattice version of the homoge- 
neous magnetic field. This has already been done before in |L4| in connection 



with the abelian Higgs model, or (2+1) QED |Tj|. We follow a slightly different 
prescription, which we describe below. 

Since we would like to impose an external homogeneous magnetic field in 
the z direction, we choose the external gauge potential in such a way that the 
plaquettes in the xy plane equal B, while all other plaquettes equal zero. One 
way in which this can be achieved is through the choice: A z (x,y,z) = 0, for all 
integers x,y,z, running from 1 to N and 

A x (x,y,z) = ~(y -l),x? N, A x (N,y,z) = ~(N + l)(y - 1), (18) 

Ay(x, y,z) = +^(x-l),y^ N, A y (x, N, z) = +| (N + l)(x - 1). (19) 

where iV 3 is the number of points on the (cubic) lattice and the coordinates 
x,y,z are integers running from 1 to N. It is trivial to check out that all pla- 
quettes starting at (x, y,z), with the exception of the one starting at (N,N,z), 
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Figure 7: Compact and non-compact L t versus ( for (3 gs = 6.0 and (3hs = 0.700. 

equal B. The latter plaquette equals (1 — N 2 )B = B — (N 2 B). One may say 
that the flux is homogeneous over the entire xy cross section of the lattice and 
equals B. The additional flux of —(N 2 B) can be understood by the fact that the 
lattice is a torus, that is a closed surface, and the Maxwell equation V • B = 
implies that the magnetic flux through the lattice should vanish. This means 
that, if periodic boundary conditions are used for the gauge field, the total flux 
of any configuration should be zero, so the (positive, say) flux B, penetrating 
the majority of the plaquettes, will be accompanied by a compensating negative 
flux —(N 2 B) in a single plaquette. This compensating flux should be "invisible", 
that is it should have no observable physical effects. This is the case if the flux 
is an integer multiple of 2ir : N 2 B = n2n — > B — nj^, where n is an integer. 
Thus we may say (disregarding the "invisible" flux) that the magnetic field is 
homogeneous over the entire cross section of the lattice. The integer n may be 
chosen to lie in the interval [0, ^j-], with the understanding that the model with 
integers m between ^ and N 2 is equivalent to the model with integers taking on 
the values N 2 — n, which are among the ones that have already been considered. 
It follows that the magnetic field strength B in lattice units lies between and tt. 
The physical magnetic field B p h ys is related to B through B = g s a 2 B p h ys , and the 
physical field may go to infinity letting the lattice spacing a s go to zero, while B 
is kept constant. 

One of the first questions which may be asked in connection with the existence 
of the external filed is of course its penetration in the bulk of the lattice. It is this 
question which is treated in table 2. The first column in the table is the integer 
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determining the magnetic field; according to the previous discussion, its range is 
from to 32 for the 8 3 lattice that we have used (in table 2 we show the results 
up to n = 20 only). Apart from P s and Pt we show the expectation value P S (B) 
of the space-like plaquette which includes the external field. Its definition reads: 
P S (P>) =< to J2x cos(F 12 (x) + B) >. The last column contains the number Q 
of vortices that penetrate into the lattice to shield the external magnetic field 
(Meissner effect). For this topological number Q we use the definition of ||14j| . We 
consider the space-like plaquette P s (x) and define the quantity P s (x), and the 
topological number q s (x) for the space-like plaquette at position x through the 
equality: 

P s (x) = P s (x) + 2nq s (x) 

and the demand that P s lies in the interval (— 7T, +ir]. The topological number Q is 
the sum of the quantities q s (x) over space-like plaquettes with fixed z-coordinates: 

xy 

We show the results for a symmetric (£ = 1.0) 8 3 lattice, with three different 
values for (3 ha . 

For j3 hs = 0.200 (which corresponds to the Coulomb phase) we observe that 
neither the space-like nor the time-like plaquettes depend on the magnetic field. 
(The two kinds of plaquettes could in principle be different even for a symmetric 
lattice, because of the presence of the external field). The same is true for < p 2 >, 
which has a constant value corresponding to the Coulomb phase. The quantity 
P S {B>) varies with the magnetic field, but its variation is trivial, in the sense that 
its value is the product of P s and cos B. We recall that: 

P S (B) = (cosB) < cosFi2 > -(sin B) < sinF 12 > . 

This suggests that the magnetic field penetrates completely and there is no sign 
of shielding. The behaviour of the topological number Q is consistent with this 
picture: it is zero for all the values of the magnetic field, indicating that no vortices 
are formed by the fluctuating gauge field and the background field penetrates with 
no obstacles into the bulk of the lattice. 

For (3hs = 0.500 the system is in the Higgs phase, as may be seen from the 
value of < p 2 > . For relatively small values of the external field (up to n = 2) P t 
and P S (B), as well as < p 2 >, are approximately constant, while P s changes even 
for small fields. It seems that the system tunes F12 in such a way as to compensate 
the change of B, so that P S (B) is constant. Small magnetic fields appear to get 
shielded by the dynamics of the scalar and gauge bosons. The picture changes 
for bigger external fields: P s changes and moves towards a value characterizing 
the Coulomb phase and this phenomenon is also confirmed from < p 2 > . The 
picture is that too big magnetic fields cannot be shielded any more and they 
penetrate into the lattice; then we pass to the Coulomb phase without Meissner 
effect. When we pass to the Coulomb phase P S {B) satisfies P S (B) = (cos B)P S 
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approximately. One may check n — 6. This picture is further supported by the 
variation of Q. This number is 1 for n — 1„ which signals the creation of a vortex 
to exactly cancel this magnetic field, while for n = 2 it equals 2, shielding the 
external field completely. For n = 3 and n = 4 no new vortices are created, so 
the screening is incomplete, and finally, for even bigger n, even these two vortices 
disappear. 



Table 2 




Volume = 


= 8 3 , C = 


= 1.0, Pgs - 


= 6.0, x 


= 2.0 




n 


cosB 


Ps 


P S (B) 


Pt 


P 2 


Q 


Phs = 0.200 





1.0 


0.9429 


0.9248 


0.9429 


1.03 





4 


0.9238 


0.9428 


0.8711 


0.9429 


1.03 





6 


0.8314 


0.9427 


0.7838 


0.9429 


1.04 





8 


0.7071 


0.9429 


0.6667 


0.9429 


1.03 





12 


0.3826 


0.9427 


0.3608 


0.9429 


1.04 





16 


0.00 


0.9427 


0.0000 


0.9429 


1.04 





20 


-0.382 


0.9428 


-0.3608 


0.9429 


1.03 





Phs = 0.500 





1.0 


0.9526 


0.9526 


0.9527 


12.08 





1 


0.9951 


0.9479 


0.9525 


0.9527 


12.09 


1 


2 


0.9807 


0.9338 


0.9520 


0.9526 


12.09 


2 


3 


0.9569 


0.9226 


0.9367 


0.9518 


11.40 


2 


4 


0.9238 


0.9250 


0.9003 


0.9506 


10.31 


2 


5 


0.8819 


0.9325 


0.8235 


0.9493 


9.05 





6 


0.8314 


0.9325 


0.7765 


0.9488 


8.54 





8 


0.7071 


0.9332 


0.6612 


0.9480 


7.72 





12 


0.3826 


0.9346 


0.3585 


0.9464 


6.01 





16 


0.00 


0.9376 


0.0004 


0.9455 


4.82 





20 


-0.3826 


0.9394 


-0.3594 


0.9448 


3.90 





Phs = 0.700 





1.0 


0.9570 


0.9570 


0.9571 


13.92 





2 


0.9807 


0.9380 


0.9563 


0.9569 


13.92 


2 


4 


0.9238 


0.8825 


0.9550 


0.9568 


13.90 


4 


6 


0.8314 


0.7920 


0.9522 


0.9562 


13.89 


6 


7 


0.7730 


0.7349 


0.9502 


0.9559 


13.90 


7 


8 


0.7071 


0.6706 


0.9475 


0.9555 


13.90 


8 


9 


0.6343 


0.6904 


0.8905 


0.9546 


13.19 


7 


10 


0.5555 


0.8011 


0.7411 


0.9535 


11.93 


4 


12 


0.3826 


0.8149 


0.5614 


0.9525 


11.17 


3 


14 


0.1950 


0.8499 


0.2684 


0.9517 


10.36 


1 


16 


0.000 


0.8561 


0.0167 


0.9512 


9.89 





18 


-0.1950 


0.8593 


-0.1531 


0.9509 


9.59 





20 


-0.3826 


0.8638 


-0.3175 


0.9507 


9.41 
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For (3hs = 0.700 the system is even deeper than before in the Higgs phase, as 
is obvious from the last column. This enables the system to screen the external 
field more efficiently: P S (B) is almost constant up to n = 8, as compared to 
n = 2 previously. The variation of P s starts early again. The penetration of B 
and the transition to the Coulomb phase is slower than before. One may check 
easily that P S (B) is no longer the product of P s and cosB, at least for small B; 
for large B, which is the region where the system approaches the Coulomb phase, 
the equation P S (B) = (cos B)P S is again approximately satisfied (one may check, 
for example, n = 16). The behaviour of the topological number is interesting 
here: the system creates Q vortices, with Q = n, for n < 8. As the magnetic 
field increases, the vortex number Q decreases and the shielding is incomplete; 
for n > 16 no vortex is present any more and the system moves towards complete 
penetration. 

Finally, we note that P t is not very sensitive to the presence of the magnetic 
field anyway. 



Table 3 






Qs(x) 


at 


random 
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To follow more closely the vortex formation, we have concentrated on n = 2 
and studied the topological numbers q s (x) characterizing each plaquette in the xy 
planes for both the S and the A models for f3 gs = 6.0, fihs = 0.700. The parameters 
are chosen such that the system is deeply in the Higgs phase; thus, we expect 
that the system will form its own vortices to shield the external magnetic field. 
These vortices will reflect the winding number of the background field, namely 
n = 2. In table 3 we show these numbers for two xy cross sections of the lattice 
(one for ( = 1.0 and one for ( = 0.1); the z-coordinates of the planes are chosen 
at random. It is obvious that the system adjusts itself, so that it screens the 
(small) magnetic field. Significantly more vortices appear in the A model; this is 
presumably due to energetics, since in the A model each layer acts on its own. 
The most important difference is that the (two) vortices appearing in the S model 
are at the same positions for all the xy planes; on the contrary, the planes in the 
A model appear decoupled and the vortices on each one of them have no relation 
to the ones of its nearest neighbours. In addition, we point out that in the A 
model the total sum of the topological numbers over the cross section need not 
necessarily yield two, since there are big fluctuations; for instance, in the case we 
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Figure 8: Time history of the topological number: ( = 1.0, (3 gs = 6.0, flhs 
0.7, n = 2. 



show in table 3, it is +1. We have also used another definition |J for the winding 
number, in which one first brings the gauge invariant link variables in the interval 
(— 7r, +7r] and afterwards follows the same steps as above; this definition has the 
advantage of being additive, so the sum over a closed surface, such as the cross 
section of a lattice should be zero. The results (not shown) have been that in the 
S model the two vortices shown before are still present, but two more vortices 
with Q = — 1 appear, so that the sum vanishes; this is also the case for the A 
model. 

We illustrate in figure 8 the evolution of the winding number Q in the S model 
(C = 1.0). We see that for the symmetric model the topological number changes 
abruptly and quickly from to n = 2. 

The picture is totally different in the A model. In figures 9 and 10 we show 
the evolution of Q for two neighbouring layers. We observe that the topological 
number changes very slowly and the final value n = 2 is reached only after passing 
slowly through intermediate values with significant fluctuations. Moreover, it is 
by no means sure that the system has settled in its final topological number and 
that no new fluctuations are going to occur. This behaviour is presumably due to 
the decoupling of the xy layers. We have found that the decoupling takes place 
at £ ~ 0.6, as one decreases ( from one to zero. 

In figures 11 and 12 we have shown the evolution of < p 2 > for the S and the 
A model respectively. They reflect the features of the corresponding changes in 
Q shown in the figures 8, 9 and 10: a sharp and quick increase for the S model 
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Figure 9: Time history of the topological number: ( = 0.1, (3 gs = 6.0, (3h s 
0.7, n = 2. 
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Figure 10: Same as in figure 9; neighbouring layer. 
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Figure 11: Time history of p 2 for ( = 1.0, f3 gs = 6.0, f3h s = 0.7, n = 2. 

and a gradual one for the A model, since there the layers decouple and do not 
move together to the new state characterized by a non-zero topological number. 
This may be an indication for the generation of the so-called pancake vortices in 
the layered phase. 
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